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ABSTRACT 


The payload of several future unmanned space missions will include a sophisticated 
gamma ray spectrometer. Severely constrained data rates during certain portions of these 
missions could limit, the possible science return from this instrument. This report investi- 
gates the application of universal noiseless coding techniques to represent gamma ray 
spectrometer data more efficiently without any loss in data integrity. Performance results 
demonstrate compression factors from 2.5:1 to 20:1 in comparison to a standard rep- 
resentation. Feasibility was also demonstrated by implementing a microprocessor 
breadboard coder/decoder using an Intel 8086 processor. 
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NOISELESS CODING 

FOR THE GAMMA RAY SPECTROMETER 
' INTRODUCTION 

Both the Mars Observer and the Mariner Mark II Comet Rendezvous Asteroid 
Flyby (CRAF) missions are expected to fly "spectrum" instruments. Such Instruments 
must periodically communicate spectrum samples in the form of histograms. This report 
investigates the application of universal noiseless coding techniques to represent such 
data more efficiently without any loss in data integrity. While investigation focuses on 
the demanding task of efficiently coding 8192 bin Gamrna Ray Spectrometer spectra 
for the Mars Observer mission, the techniques and analysis provided should more generally 
apply to other spectrum instruments with little modification. 

SUMMARY 

A single coding algorithm which will efficiently represent gamma ray spectrometer 
spectra at any spectrum update rate from 5 seconds up to 5 minutes is identified. Data 
rate compression factors from 20:1 to 2.5:1 compared to a standard spectrum represen- 
tation were demonstrated over this range using representative Mars Observer Gamma 
Ray Spectrometer test data. Looking at these results in another way, spectrum updates 
can be accomplished 7 to 23 times more frequently with noiseless coding 
than with a standard representation if both approaches use the same data rate, R, 
(where 250 R <, 800 bits/s). Whereas spectrum updates every 30 seconds would 
require roughly 2500 bits/s by standard means, spectrum updates using noiseless cod- 
ing could be accomplished every 10 seconds using only 500 bits/s. 

The presentation of these results will begin with discussions leading to the standard 
data representation for instrument generated spectra to be used as a means of comparison. 
We will then develop the structure and performance criteria for a Gamma Ray Spectrometer 
Noiseless Coder and finally graph and tabulate the results of performance runs on a 
representative test set. Additional techniques are developed to improve performance under 
unusual burst mode conditions. 


1 



THE BASIC INSTRUMENT 


A "spectrum instrument" monitors the rate of arrival of photons or particles which 
have energies lying within one or more instrument energy ranges. This is accomplished 
by counting the number of hits in each "energy bin" over a known sampling interval. 
The average rate for a given energy bin is simply the number of hits divided by the sam- 
pling interval. A spectrum can be viewed as the composite of the counts or arrival rates 
of all the instrument energy bins (essentially a histogram). 

"Raw" Data Representation 

Let 


N 0) 

be the number of possible energy bins detectable by the spectrum instrument. And |;;tt 

n = T log2N ~] bits (2) 

be the number of bits required in a fixed length binary code to identify any one of the 
N bins. Of practical significance is the case where 

N = 2 n 


as for the proposed Mars Observer Gamma Ray Spectrometer Instrument where n = 1 3 
bits and N = 8192 bins. 

Now assume that all the hits occurring in a time interval of 

r seconds (3) 

are collected in a buffer in order of occurrence. Letting 

a (4) 


t f x”| is the smallest integer greater than or equal to x. 
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be the number of such hits and a; the bin number of the I™ hit, the buffer contents 
(prefaced by the value of ct T ) would appear as the sequence 

B = a T 8i a 2 a 3 ••• a a • (5) 

T 

S' is viewed as the “raw data" mode of representation for a spectrum instrument. By 
communicating B every r seconds we will know the order of occurrence of each hit in 
the interval but we will not know precisely when they occurred. 

Consider the raw data rate. Note that the number of hits in an interval of length r 
is not deterministic. The number of arrivals is in fact well modeled as Poisson^] with 
parameter 


X hits/s (6) 

as the composite average arrival rate for all the energy bins. Then, the expected number 
of hits in interval r is given as 

E{« t ) = Xr hits . (7) 


For a given Xr, we use 


L bits (3) 

to represent a r where any number greater than 2*- - 1 is set to 2^- - 1 and the remaining 
samples of a T are discarded. L is chosen so that there is little chance that 2*-- 1 is 
exceeded. 

The total average bit rate required to communicate B in (5) is then 



nXr + L 
r 


= nX + L/r bits/s 


(9) 
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where we see that the contribution for L/r becomes negligible os r increases. 

For the Mars Observer primary operating environment wo hove X = 250 hlts/s, 
n = 13. If we take t <=< .05 seconds, the expected number of hits is 12.5. By using 
L = 8, the number of hits would hove to exceed 255 before any loss In data occurred. 
The required "raw" data rate is then, from (9), 

R r b 3410' bits/s . (10) 

Under rare burst mode conditions (solar flares) X may Increase by an order of magni- 
tude, Then taking L = 1 2, we have 

R R = 32740 bits/s (11) 


under burst mode conditions. 
Standard Spectrum Representation 
Let 


S — 8^2 Sjyj (12) 

be the complete sampled spectrum generated over a t second Interval where 

I count of the 

number of hits (13) 

in energy bin i 

We will later focus on minimizing the average number of bits needed to represent S for 
any choice of t. Here we will define the standard representation of S which will later 
be used as a means of comparison. 

Let mj(r) be the number of hits in bin i in a time interval of r seconds. If by observa- 
tion or probability model we can say that 

Pr [m j (r) > M, max (r)] < e (14) 

where e is very small, we need no more than 
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bj<r) = riog 2 M l maX ( T n bits 'im 

to represent almost any occurrence of hits in bin i in an interval of r seconds. As a 
backup, if a count were larger than 2^^ - 1 (with probability < e) we could, at worst, 
truncate the Identification of number of counts to 2 b i tT '- 1. 

Obviously, a smaller interval, r, and or a lower hit rate means there Is less chance 
of the number of hits exceeding some threshold. Thus to achieve a given level of protec- 
tion from count truncation (e in (14)) the number of bits required in (15) can be quite 
different for different I and t. This is particularly true for the Gamma Ray Spectrometer. 
However, standard spectrum Instrument data handling does not capitalize on the broad 
disparity in quantization requirements across all the bins and possible spectrum inter- 
vals. Basically, the quantization requirements for all bins, b*, is chosen as 

b* = M £X pog^H^lrn . (16) 

That is, b* satisfies the worst case requirements for ail bins and the longest expected 

spectrum interval. The communication of each spectrum then requires 

b* bits/bin . (17) 

By (1) there are N bins so that spectra communicated every t seconds require 

Nb* 

r ( T ) = — bits/s (18) 

o r 

which we define as the “Standard Spectrum Data Rate". 

A plot of Rg(r) and Hr is shown in Fig. 1 where we have taken 

73700 

b* = 9 so that Rq(t) = — - — bits/s (19) 

for the Mars Observer Gamma Ray Spectrometer. 

We have marked some particular points of interest on the graph: 

1 ) The standard spectrum representation is better than a raw data representation 
at spectrum update intervals in excess of about 21 seconds. 
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Fig. 1. Data Rates for Standard Spectrum Representation, Rg(r), 


2) Spectrum updates as frequently as every 30 seconds are highly desirable from 
a science point of view since update rate determines ground resolution on Mars 
Observer. Spectrum updates at 30 second intervals would require almost 2500 
bits/s which is an order of magnitude more than may be available during certain 
portions of the Mars Observer mission. 


STATISTICAL CHARACTERISTICS (Gamma Ray Spectrometer) 11 ] 

The arrival of hits in the i^ bin is made up of two independent components: 
background (mostly cosmic rays), and the spectral emissions from the elements and com- 
pounds which the instrument hopes to detect. The latter category we label as the science 
component. 

The arrival of hits from both components can be closely modeled as Poisson with 
parameters 
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< 20 ) 


X. B t and Xpr 

respectively. 

The probability of k hits in time interval r by a Poisson process with arrival rate X is 

P k (r) = e‘ Xr , k^0,r^0 (21) 

and the expected number of hits is 

E(k) = Xr hits . (22) 

The overall arrival process for the i 1 ^ bin is Poisson with parameter XjT given as 





(23) 


Can we consider the arrival process independent from bin to bin? The background 
process is clearly independent from bin to bin. On the other hand, a single compound 
may generate hits in several bins (characteristic lines) so that this is not entirely true for 
the science component. However 


X. B > > Xp (24) 

and the coding process we will consider later will deal with intervals much too short to 
distinguish spectral lines (if they were even there; the science component may or may 
no; be present). Then for all practical purposes we can take the overall arrival process 
as Poisson and independent from bin to bin. 

Note that by the same arguments, the arrival of instrument hits is Poisson with 
arrival rate 


N 


X = 



where X is the same term introduced in (6). 


( 25 ) 
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Observations on the \| 

The following observations will influence our spectrum coding approach in subse- 
quent sections. 

• The Xf's tend to decrease with increasing bin value i. (26) 

• The largest \j may be as much as three orders of magnitude larger than 

the smallest. (27) 

• Adjacent Xj's tend to be similar. That is, if i =* j then Xj = Xj . (28) 

We are now ready to develop an efficient spectrum noiseless coder. 
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II. SPECTRUM NOISELESS CODER 


NOTATION AND DEFINITION 

The application of noiseless code operator 

(29) 

to sequence X produces the coded sequence t^X], and given MX] we can retrieve the 
original X precisely. That is, we can recover the original "noiselessly". 

In general we subscript and superscript \f/[»] to specify different choices of operators. 
The particular operator then means those operations that are performed to produce a 
coded sequence. 


£(MXj) (30) 

means the "length" of sequence \f/[X] so £ [•] means the calculation of the length of some 
sequence. 

When emphasizing the concatenation of two or more items we will use an asterisk, *. 


PREPROCESSING 

The basic form for a noiseless coder is shown in Fig. [3], 

The first steps are to perform instrument-specific "preprocessing" operations which 
seek to produce a sequence of independent non-negative integers, 5, such that a smaller 
integer tends to occur more often than a larger one. That is, 

Pr[S = 0] 5: Pr[5 = 1] 5; Pr[5 = 2] £ . . . . (31) 

The second instrument-independent step, indicated by ^ [•] in Fig. 2, is to efficiently 
represent such sequences of integers. The underlying technique is to assign shorter code 
words to the more likely smaller integers and larger code words to the less likely larger 
integers such that, on the average, fewer bits are used. We will not need to worry about 
designing ^ [•]. References 2 and 3 document a broad class of adaptive algorithms (code 
operators) which will provide efficient coding of a preprocessed 8 sequence under vir- 
tually any practical situation. We will later merely identify the operator we elect to use. 
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I J 


Fig. 2. Basic Noiseless Coder Structure. 

The problem here is to identify the instrument-specific preprocessing that is appro- 
priate for the Gamma Ray Spectrometer. 

Prediction and Ordering 

Letting xj be the prediction of the actual number of hits Xj, the prediction error is 
given as 

A. = ^ - x, . (32) 

Then assuming we know Xj for the i^ bin we take x; as the closest integer to the 
expected number of hits, that is from (22) t and (23) 

^ = [Xjr] . (33) 

The resulting error distribution is given directly by the Poisson distribution 

t[a] is the closest integer to a. 
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Pr[A, = j] * Pr[k = [X.t] + j] for any j 


(34) 


where the latter term can be calculated by substitution In (21). Using this distribution 
we can got the exact ordering which would allow us to map the error Aj into the Integers 
such that condition (31) is met. 

However, a prediction about the mean, Xj-r, generally results in a unimodal error 
distribution around zero. When XjT >> 1 the standard error mapping in Table 1 gener- 
ates the same ordering (for ail terms of consequence) as using the Poisson distribution 
itself. This can be modified to account for the boundary conditions (when XjT is near zero) 
resulting from the fact that the number of counts cannot be less than zero^J. 


Table 1. Basic Mapping of Aj Into the Integers 5|. 


PREDICTION 

ERROR 

A i 

INTEGER 

*1 

0 

0 

+ 1 

1 

-1 

2 

+ 2 

3 

-2 

4 

+ 3 

5 

-3 

6 

• 

• 

• 


Specifically, with x^^x equal to the largest allowed count value, we can map the 
error Aj as 
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2A, - 1 


if 0 < A, ^ S) 
if A| > x. 

( 36 ) 

if x r x MAX s A i ^ 0 
lfx i“" x MAX >A i • 

This mapping yields almost identical results, under all conditions, as that produced 
by direct use of the Poisson distribution. 

Based on these observations, all future discussions will presume one of the follow- 
ing preprocessing options: 

* We do not attempt to predict at all and simply take 

5. = Xj (36) 

for each sample (bin). 

• We predict Xj as the closest integer to Xjr and then use the mapping in 
(35) to convert the prediction errors into integers. The Xjr are derived from 

long term average statistics. (37) 

PERFORMANCE MEASURE, ENTROPY 

Consider a single bin spectrum driven by a Poisson process with parameter Xr. Then the 
probability of k hits, p k (r), is given by Eq. 21. The entropy of this distribution is given as 

00 

H(Xr) = - P k (r) log 2 P k (r) . (38) 

k -0 

This expression represents the minimum possible average rate (averaged over many 
spectra) that any noiseless coder could code the counts of this single bin spectrum. A 
plot of H(Xr) vs. Xr is given in Fig. 3. 


214,1 

x MAX " x i 
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ENTROPY. BITS/ BIN 



Fig. 3. Poisson Entropy, 


If \t were fixed and known then a long sequence of counts taken at r second inter- 
vals could be represented with no fewer than H(Xr) bits/bin count, A coder which coded 
closely to H(Xr) would be considered efficient. 


Now suppose we applied a separate noiseless coder to each of the N = 8192 inde- 
pendent bins of the Gamma Ray Spectrometer. We could expect to code no better than 


H x !rl = 


N 

Y] H(Xjr) bits/bin 
i = 1 


(39) 


Using long term average values for the (Xj) we call this term the Gamma Ray Average 
Entropy. 


On the other hand, if we restricted ourselves to a single fixed noiseless coder for all 
bins, we could expect to code no better than 


H 



£ H x (rl 


where X is the average arrival rate over all bins given in (24). 


( 40 ) 
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For the Gpmma Ray Spectrometer we must also worry about efficiently coding each 
of the different bins for a broad range of r. The performance goal we're shooting for is 
to code close to H^(t) in (39) for any r of interest. 

A plot of the various entropies derived from representative Mars Observer test data 
Is shown in Fig. 4. 


These same curves have been converted to the equivalent data rate by the formula 
(see 18) 


data rate = 


rate in' 
bits/bln 


bits/s 


T 


(41) 


and are plotted again in Fig. 5. H(Xr/N) and H^(t) have been replaced by H^Xr/N) and 
H^(r) respectively. 



Fig. 4. Spectrum Coder Entropy Performance Bounds. 
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DATA RATE. BITS/s 



Fig. 5. Spectrum Coder Performance Bounds in bits/s. 

The first observation from these curves is that the lower curve for H^(t) is dramati- 
cally lower than the corresponding rate for a standard representation. Second, the curve 
for Hx(r) is significantly lower than that of H'(Xt/N). Unfortunately, the apparent implica- 
tions for achieving the performance of the lower curve is to build a separate coder for 
each \t combination. This will not be necessary. In the following sections we will com- 
plete the definition of a coder which approaches the performance of the lower curve, 
for all 7 , with complexity closer to that of a single coder. 


COMPLETING THE DEFINITION 
Partitioning the Spectrum 

The observations in (25)-(27) give us the necessary clues to avoid the implied indi- 
vidual coder for each bin (and r). We need only partition the spectrum into subblocks and 
code each subblock separately. So let S = s-jS 2 S 3 . . . . S|\j in (12) be rewritten as 

S = X 1 *X 2 *X 3 * . . . *X M (42) 
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where the Xj's are consecutive subbiocks of the samples of §. That is, if s v is the last sam- 
ple of X| then s v + <| is the first sample of the next subblock, Xj + -|. 

Gamma Ray Code Operator 

We define the overall Gamma Ray Spectrometer spectrum noiseless code operator 
as i Pq[»] where 

* G (S] = ■ y>W (43) 

and where ^ g [«] is a code operator that works on subblocks. In this report we take 
as the structure in Fig. 2 with the preprocessing assumptions in (36) and (37) and 

^[«] = ^ 10 H (44a) 


defined explicitly in Ref. 2 or 

M = 8-option FAST Compressor (44b) 

defined in Ref. 3. Roth of these are too involved to go into here. 

The next section will look at the performance of this configuration. 

Breadboard 

A breadboard coder/decoder of (using = i/'ioM) was recently built, and 
tested by Bruce Parham using an Intel 8086 microprocessor,^] The breadboard coder 
requires roughly 1 Kbyte of memory for instructions and another 1 Kbyte for internal buf- 
fers. A decoder is about half this size. Both will handle throughput data rates of 20 
kbits/s, far higher than the allowable data rates on Mars Observer. 


SPECTRUM CODER PERFORMANCE RESULTS 
Using ^ 10 M 

Gamma Ray Spectrometer test data, highly representative of the expected Mars 
Observer environment, was provided by Jack Trombka of Goddard^], it consisted of 
roughly 1/2 hour of raw data acquired from a recent balloon flight. 
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The noiseless coder just defined was run on the entire data set for different values 
of r, running from 6 seconds up to 6 minutes 1 '. The performance results are displayed 
in Figs. 6-9. 

Figure 6 shows performance in bits/bin with the Gamma Ray Spectrometer noiseless 
coder rate denoted by Rq(t) and where 

Rq(t) « 0.189 Vr bits/bin . (461 



Fig. 6. '/'qM Performance in bits/bin Using ^qM. 


tpach of these curves presumes the preprocessing assumption in (37). Using the simple 
approach in (36) (basically no a priori knowledge) resulted in only minor losses. Perfor- 
mance reductions were only observed at data rates of less than 500 bits/s and were 
typically less than 4%. The largest performance loss observed was less than 10%. 
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Fig, 7. 0 q[*] Performance in bits/s Using ^qH’ 
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Fig. 8. Data Rate Compression Factor, CF r (t) Using 
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Fig, 9. ^ Q [»] Spectrum Update (r) Compression Factor, CF r (R) Using l^gM. 


These curves have been translated to bits/s in Fig. 7 using the relation in (41). The 
corresponding performance for the Gamma Ray Spectrometer noiseless coder is denoted 
by R&(r) where 


Rq<t) « - 1 bits/s . 

° VT 


(46) 


Additionally, the performance of relative to a standard representation is 

shown in Fig. 8 as a data rate compression factor CF R (r) given by 


CF R (t) 


9 

R g <t) 


R s (r) 

RqW 


47.6 


(47) 


when \^oM * s usec *- 
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Figure 8 compares the average communication rates that can be achieved for a 
given spectrum update rate. Perhaps more important is a comparison of the spectrum 
update rates that can be achieved for a given data rate. For this we define a second com- 
pression factor CF r (R), where 

r such that the 
standard spectrum 
representation 
requires R bits/s 

(48) 

t' such that ) 

performance > 

requires R bits/s J 

which is very close to a straight line with intercept at R = 0 when ^iqM is used. Under 
these assumptions we have 

ciym » ■ (49i 

In fact it is this observation that led to the simpler but accurate approximations in 
(45) - (47). CF r (R) is plotted in Fig. 9. 

Observations. Although these figures need to be studied closely, we can make 
several important observations. 

• The ^q[»] coder is efficient. Its performance is close to the optimistic perfor- 
mance bound for multiple coders, H^(r). 

• Data rate compression factors at a given spectrum update interval, r, range from 
20:1 to 2.5:1 for 5s <, r <, 300s. 

• Compression factors for spectrum update intervals at a given data rate, R, range 
from 7.3:1 to 23:1 when data is constrained to 250 < R <, 800 bits/s. 
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Of particular note are the following: 

• Whereas spectrum updates every 30 seconds would require roughly 2500 bits/s 
by standard means, spectrum updates using noiseless coding could be ac- 
complished every 10 seconds using only 500 bits/s. 

• A data rate of 500 bits/s would allow spectrum updates no more often than once 
every 1 50 seconds if a standard representation were used. This Is an update 
interval compression factor, CF r (R), of 15:1. 

• i A q[ * ] c °dlng data rate, Rq(t), with an update interval of 1 50 seconds drops to 
100 bits/s. 

• v^gM can code spectrums every 30 seconds using only 300 bits/s compared 
to approximately 2500 bits/s by standard means. This is a data rate compres- 
sion of CFp(r) = 8.3:1 . 

• A 300 bits/s data rate would allow spectrum updates no more frequently than 
once every 230 seconds. This is a spectrum update interval compression of 
CF r (R) = 7.6:1. 

Using a Simple Operator 

Code Operator ^iqM includes special modes to deal with data characteristics 
corresponding to very low entropies. For the Gamma Ray Spectrometer spectrums this 
situation occurs at the top end of a spectrum (small X) when r is small (less time to ac- 
cumulate hits). A simpler code operator is feasible if one is willing to give up some per- 
formance at low values of r. 

An 8-option "FAST” compressor^ 3 ] operating on blocks of 10 bins was applied 
to the same test set. t We will denote this simpler algorithm by 

. (50) 

A performance comparison of i^p[*] with ^iqM is shown in Figs. 10-12. As the 
figures show, the advantage of i/'iq(«] becomes significant at spectrum intervals less 
than 1 minute. At r = 30 s, ^p[»J requires about 30% more data rate than ^iq[»). The 
advantage of \£"|qM increases rapidly with smaller r. 


^This is the same algorithm which will be used to compress images at the 1986 
Voyager II encounter of Uranus. 
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Fig. 12. i Pq[»] Data Rate Compression Factor, CF R (r), Using ^qM or 



III. BURST MODE CODER 


Solar flares may cause a brief increase in the average arrival rate by as much as 
an order of magnitude. Under these conditions it is desirable to sample spectrums at up- 
date intervals as low as 0.05 s/spectrum. 

The standard spectrum representation would require over 1 Mbit/s to accommodate 
this situation. On the other hand, the "raw data" mode in (5) is far more efficient, re- 
quiring, from (1 1 ) 


R r = 32740 bits/s 


(51) 


ESTIMATING ^ G [*] PERFORMANCE 

We can us8 previous results to estimate the performance of i/'qM under these 
burst conditions. 

Let X N and X^ , i = 1,2, ..., 8192 represent arrival rates under nominal conditions 
and X* 3 and Xj b , i = 1,2, ..., 8192 represent arrival rates under burst mode condi- 
tions. Further, let r b denote the burst mode spectrum sampling interval. 

Now if we satisfy 



(52) 


by assuming 


xf = 10 X^ for all i 


(53) 


we would expect the same "bits/bin" performance from X G [*J when 

| X = X^ , r = as when |x«=X^,t=10t^| 

because, in both cases, each bin should (statistically) receive the same number of hits. 
From Eq. 6 the expected rate under burst mode conditions becomes 


24 



Rq (r b ) = 0.189 \/l0 7 b bits/bin 


(54) 


Using Eq. 41 this translates to data rate as 




18192)10.189) /To 


4896 


bits/s 


(55) 


Substituting, we find that i^ G (*] spectrum coding requires 

R b ' (0.05) = 21896 bits/s (56) 

under the stated burst mode conditions. This is roughly 1 1 ,000 bits/s less than required 
by a "raw data" representation (51). Subsequent paragraphs wii! seek to improve this 
result. 

CODING THE DATA DIRECTLY 

Rather than first collecting a spectrum we will consider coding the row data sequence 
in (5), repeated here for convenience. 

s = “,‘ 1 3 2 ’"'\ 1571 

where a T is the number of hits in r seconds and aj the value of the i™ hit. 

By earlier discussions the contribution of a T to total rate is negligible compared to 
the sequence 

A — a,j 3^ .... a a (58) 

so we will concentrate on A. 

Since we are ultimately interested in a spectrum from A, the order of arrivals is 
unimportant. Then we define the operation SORT[«] as a sorting algorithm such that 

SORTfA] = C = c 1 c 2 . . . . c (59) 
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where the Cj are the samples of A, rearranged to place any aj in front of any larger aj. 
That is, 

Cj C| if i < j . (60) 

Now form the sequence 

d = diff;cj = d 1 d 2 . . . . d 

where 

d 1 = c 1 ^ 0 

and (61) 

dj = Cj - C|__ 1 ^ 0 for i > 1 . 

We now define the coder of C by 

^[C] = a r *^D] (62) 

where it is only necessary to select a proper coder for D. 

We investigated the following noiseless coder for D 

\M‘l = = FAST COMPRESSOR (see (50) and Ref. 3): 

• operating on blocks of 10 samples ; 

• using 13 standard fast compressor modes; and (63) 

• one special mode which identifies the presence of an all-zeroes 
block. 

A block diagram of the overall coder of B is shown in Fig. 13. 

It should be understood that the SORT process is not reversible. Only C can be re- 
covered precisely. However, a spectrum can be derived without error from C. 
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Fig. 13. Coder of Raw Data Sequence. 


Performance 

For notational purposes, let 


r'(r b ) (64) 

denote the data rate required by the burst mode coding algorithm i^ r [ •], in (59)— (62) and 
Fig. 13, when the interval between the transmission of collected hits is once every 
r = r b seconds. A performance comparison of the burst mode coder and the spectrum 
coder is shown in Fig. 14 (under burst mode conditions). 

Performance of the burst mode coder is better when < 0.25 s. At the desired 
interval of = 0.05 s, ^ r [«] requires 

r'(r b = 0,05) =16 kbits/s (65) 

compared to about 22 kbits/s for This is a 2:1 reduction compared to a direct 

transmission of raw data (51). 
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Fig. 14. Burst Mode Coding. 
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